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Medium modification of dihadron fragmentation functions due to gluon bremsstrahlung induced 
by multiple partonic scattering is studied in both deep-inelastic scattering (DIS) off large nuclei and 
high-energy heavy-ion collisions within the same framework of twist expansion. The modified frag- 
mentation functions for dihadrons are found to follow closely that of single hadrons leading to a weak 
nuclear suppression of their ratios as measured by HERMES in DIS experiments. Meanwhile, a mod- 
erate medium enhancement of the near-side correlation of two high transverse momentum hadrons 
with increasing centrality is found in heavy-ion collisions because of the trigger bias and the in- 
crease in parton energy loss with centrality. Successful comparisons between theory and experiment 
for multi-hadron observables in both confining and deconfined environments offers comprehensive 
evidence for partonic energy loss as the mechanism of jet modification in dense matter. 
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Medium modification of the jet structure has emerged 
as a new diagnostic tool for the study of partonic prop- 
erties of the dense matter The modification goes 
beyond a mere suppression of inclusive spectra of leading 
hadrons and can be extended to include the modification 
of many particle observables, the simplest of which are 
two-hadron correlations within the jet cone. Such two- 
hadron correlations have been measured both in DIS 0| 
and high-energy heavy- ion collisions 0, While the 
two-hadron correlation is found slightly suppressed in 
DIS off a nucleus versus a nucleon target, it is moderately 
enhanced in central Au + Au collisions relative to that in 
p + p. This is in sharp contrast to the observed strong 
suppression of single inclusive spectra [B|, Q in either DIS 
or central A + A collisions and constitutes a necessary 
test of any energy loss formalism used to compute the 
modification of the single inclusive spectrum. 

Theoretically, multi-particle correlations from jet frag- 
mentation and their medium modification can be studied 
through n-hadron fragmentation functions which can be 
defined as the overlapping matrices of partonic field op- 
erators and n-hadron states. These n-hadron fragmenta- 
tion functions are non-perturbative and involve long dis- 
tance processes. However, they may be factorized from 
the hard perturbative processes and their evolution with 
momentum scale may be systematically studied in per- 
turbative QCD (pQCD) 0, similarly as the usual single- 
inclusive hadron fragmentation function (SFF) Q. In 
this respect, they serve as an important test of the fac- 
torized pQCD formalism currently being applied in the 
computation of jet-like observables. 

In this Letter, we report our first study of the medium 
modification of dihadron fragmentation functions (DFF) 
in both cold nuclei and hot deconfined partonic matter. 
The medium modification of DFFs in DIS off nuclei will 
be derived within the framework of generalized factor- 



ization and twist expansion 0, [lfl | . The results are then 
extended to the case of parton propagation in heavy-ion 
collisions. With exactly the same parameters determined 

ToL 11 1, the medium modifica- 



from the modified SFFs 
tion of two-hadron correlations in both DIS off nuclei 
and heavy-ion collisions are predicted and compared to 
the experimental data. 

Similar to SFFs, the DFF of a quark into two hadrons 
of flavor hi,h 2 , with total forward momentum p~^ = p^ + 
p 2 (total momentum fraction z = z\ + z 2 ) and relative 
transverse momentum q± = pu_ — p 2 ± can be defined as 
the Fourier transform of the overlapping matrices of the 
quark fields ip q (x),ip q (§) and two-hadron inclusive final 
states \pi,p 2 ,S — 2) as, 
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and can be factorized from the hard processes 0]. They 
also satisfy the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations as derived in 
Ref. 0] • One unique feature of these equations is the con- 
tribution from independent fragmentation of two partons 
after the parton splitting. These equations have been 
solved numerically Q and the Q 2 evolution agrees very 
well with results from JETSET [l3] Monte Carlo simula- 
tions of e + + e~ — > h\ + h 2 + X processes. Although both 
the single and dihadron fragmentation functions evolve 
rapidly with Q 2 , their ratio has a very weak Q 2 depen- 
dence. In the absence of experimental data, JETSET 
Monte Carlo results will be used as the initial condition 
for the vacuum DFF in this study. For SFFs, the BKK 
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parameterization j 1 31 ] will be used; which also agrees well 
with JETSET results. 

Applying factorization to dihadron production in sin- 
gle jet events in DIS off a nucleus, e(L{) + A(p) — > 
e(L 2 ) + hi(pi) + h 2 (p2) + X , one can obtain the dihadron 
semi-inclusive cross section, 
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in terms of the semi- inclusive tensor at leading twist, 
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In the above, i A1 ^(l/2)Tr(^i 7 ^27^), the factor H» v 
represents the hard part of quark scattering with a 
virtual photon which carries a four-momentum q = 
[— Q 2 /2q~, q~ 7 0j_] and f^{x,Q 2 ) is the quark distri- 
bution in the nucleus which has a total momentum 
A[p + ,0,0^]. The hadron momentum fractions, Z\ = 
Pi/q~ and z 2 = p 2 / q~ , are defined with respect to the 
initial momentum q~ of the fragmenting quark. 

At next-to-leading twist [3, the dihadron semi- 
inclusive tensor receives contributions from multiple scat- 
tering of the struck quark off soft gluons inside the nu- 
cleus with induced gluon radiation. One can reorganize 
the total contribution (leading and next-to-leading twist) 
into a product of effective quark distribution in a nucleus, 
the hard part of photon-quark scattering H^ v and a mod- 
ified DFF: Dq uh *(zi,Z2) (which includes, within it, the 
effect of subsequent scattering and gluon radiation). The 
calculation of the modified DFF at the next-to-lcading 
twist in a nucleus proceeds [H| similarly as that for the 
modified SFFs [l(| and yields, 
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In the above, x B = -Q 2 /2p + q ,x L = l 2 _/2p + q y(l-y), 
is the transverse momentum of the radiated gluon, 
AP q —> q g and AP q ^ gq are the modified splitting functions 
with momentum fraction y, whose forms are identical to 
that in the modified SFF [k]]. The switch (hi -> h 2 ) is 
only meant for the last term, which represents indepen- 
dent fragmentation of the quark and radiated gluon. The 



corresponding modified splitting function, 
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is similar to AP q ^ qg but does not contain contributions 
from virtual corrections. In the above, Ca = N c = 3, and 
is the average intrinsic parton transverse momentum 
inside the nucleus. 

Note that both modified splitting functions depend on 
the quark-gluon correlation function T qg in the nucleus 
that also determines the modification of SFFs [13]. For 
a Gaussian nuclear distribution, it can be estimated as, 

T q A g (x B ,x L ) = C(Q 2 )m N R A f q A (x B )(l - e"^^), (6) 

where, xa = I/tunRa, ttln is the nucleon mass and 
Ra = 1.12A 1 / 3 is the nuclear radius. The small variation 
of the gluon correlation function is neglected; its average 
value is absorbed into the overall constant C. This is the 
only parameter in the modified DFF which might depend 
on the kinematics of the process but is identical to the 
parameter in the modified SFF. In the phenomcnological 
study of the SFFs in DIS off nuclei, C = 0.006 GeV 2 was 
determined within the kinematics of the HERMES ex- 
periment Q . The predicted dependence of nuclear mod- 
ification of the SFF on the momentum fraction z, initial 
quark energy v = q~ and the nuclear size R a agrees very 
well with the HERMES experimental data [ll | . 

With no additional parameters in Eq. ^ , one can pre- 
dict the nuclear modification of DFFs within the same 
kinematics. Since the DFFs are connected to SFFs via 
sum rules , it is more illustrative to study the modifica- 
tion of the distribution for the second rank hadrons nor- 
malized by the number of leading hadrons (i.e., hadrons 
with z > 0.5), 



N 2h (z 2 ) 



1-Z 2 



dz 1 D' q l ^ h2 (z 1 ,z 2 ) 



0.5 




dziD^(zi),(7) 



0.5 



where z\ and z 2 < z\ are the momentum fractions of the 
triggered (leading) and associated (secondary) hadrons, 
respectively. Shown in Fig. [T] is the predicted ratio of the 
normalized associated hadron distribution in DIS off a 
Nitrogen (A = 14) and Krypton (A = 84) target to that 
off a proton (A = 1) i.e., R 2h (z 2 ) = Nf h ( z 2 )/Nl(z 2 ), as 
compared to the HERMES experimental data [2| . 

The agreement between the prediction and the data 
is remarkable given that no free parameters are used. 
The suppression of R 2 h(z 2 ) at large Z2 with atomic num- 
ber, is quite small compared to the suppression of the 
SFFs ahj. Since N 2 h(z 2 ) is the ratio of double and sin- 
gle hadron fragmentation functions, the effect of induced 
gluon radiation or quark energy loss is mainly borne by 
the single spectra of the leading hadrons. At small val- 
ues of 22, the modified DFF rises above its vacuum coun- 
terpart more than the modified SFF. This is due to the 
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new contribution where each of the detected hadrons em- 
anates from the independent fragmentation of the quark 
and the radiated gluon. In the experiment, the measured 
v and Q 2 vary with Z\ and z 2 ; this has been incorporated 
in the calculation 
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FIG. 1: (Color online) Results on the medium modi- 
fication of the associated hadron distribution [B,2h(z2) = 
N 2h (z2) I r N2h(z2) , see Eq. 0] in nuclei versus its momentum 
fraction as compared with HERMES data @] in DIS with Ni- 
trogen (N) and Krypton (Kr) targets. 



In high-energy heavy-ion (or p+p andp+^4) collisions, 
jets are always produced in back-to-back pairs. Correla- 
tions of two high-px hadrons in azimuthal angle generally 
have two Gaussian peaks 0, 3 ■ Relative to the triggered 
hadron, away-side hadrons come from the fragmentation 
of the away-side jet and are related to SFFs. On the 
other hand, near-side hadrons come from the fragmenta- 
tion of the same jet as the triggered hadron and therefore 
are related to DFFs. 

The near-side correlation, background subtracted and 
integrated over the azimuthal angle [TBI, [l6| , can be re- 
lated to the associated hadron distribution or the ratio of 
the DFF to the SFF, both averaged over the initial jet en- 
ergy weighted with the corresponding parton production 
cross sections. We report on the first such calculation 
in this Letter. Assuming a factorization of initial and 
final state effects (as done in the case of single inclusive 
observables [ill flTl ] ) , the differential cross-section for the 
production of two high px hadrons at midrapidity from 
the collision of two nuclei A and B at an impact param- 
eter b between b rn i„, b rnax is given as, 
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where, (x a , Q 2 )[Gf (xb, Q 2 )] represents the nuclear 
parton distribution function for a parton a(b) with mo- 
mentum fractions x a (xb) in a nucleus A(B), t A (is) rep- 
resents the nuclear thickness function and da a b-^ c dj dt 
represents the hard parton cross section with Mandel- 
stam variable t. The final state momentum fractions 
z\ , Z2 represent the the momentum fractions of the two 
detected hadrons with respect to the initial jet energy. 
The factor K ~ 2 accounts for higher order corrections 
(identical to that used in the case of the single inclusive 
spectra). The medium modified DFF may be expressed 
as in Eq. , with the modified splitting functions gen- 
eralized from Eq. ([5]) as, 
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As in the case of DIS, l± represents the transverse mo- 
mentum of the radiated gluon and T M replaces the factor 
Tq g (in Eq. ([5])) and represents the parton-gluon correla- 
tion function in a quark-gluon plasma. The primary dif- 
ference between the denominators of Eqs. (O and ([9]), 
lies in the weighting of a particular jet like parton by the 
initial production process (see Ref. [l8j for the evaluation 
of T M in single inclusive observables). 

At present, there exists no measurement of high mo- 
mentum hadrons associated with a high p^ lg '(> 8GcV) 
hadron in the case of p+p collisions. Therefore, we will 
compare the associated yields in d+Au and Au+Au col- 
lisions as a function of centrality as well as the associ- 
ated pt- These yields are very sensitive to the num- 
ber of flavors detected. In Fig. the yield of charged 
hadrons with associate transverse momentum in two 
ranges (6GeV<pf? soc <p^ ig - and 4GeV<pf? soc < 6GcV), 
associated with a hadron with p^ Ig > 8GeV, as measured 
in Ref. , are presented along with three different cal- 
culations using Eqs. (HJ) and ([9]). The dashed lines cor- 
respond to the case of all charged hadrons where as the 
dot-dashed lines corresponds to the case of charged pions 
only. In both cases, the vacuum DFFs are estimated from 
JETSET Monte Carlo simulations (see Ref. 0]). Unlike 
in the case of the experimental measurements, no decay 
corrections have been introduced in these calculations: 
the number of pairs measured, includes not only those 
produced directly from the fragmentation of a jet but also 
decay products of particles from fragmentation. As a re- 
sult, such fragmentation functions tend to be somewhat 
larger at lower momentum fractions than fragmentation 
functions with such corrections. The integrated (over 
transverse momentum) associated yield being dominated 
by lower momentum fractions shows a larger effect as 
lower momentum ranges in p^ soc - are chosen. 
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A decay corrected DFF [D(zi, z 2 )] may be constructed, 
phenomcnologically, by comparing with the differential 
spectrum of associated particles in d+Au collisions (as a 
function of z T = py ig '/pf? soc -, in Ref. [3]): 



D{z u z 2 ) = 0.55 x D( Zl , z 2 )[l + lA( Zl /z 2 - 1)]. (10) 

Using the above DFF, we obtain a very good agreement 
with the data (solid lines). All estimates show a mild rise 
with centrality which originates solely from the increased 
trigger bias as the centrality of the collision is increased. 
We have focused on large p^ ls ' and p|? soc -, as large mo- 
menta ensure the validity of the independent fragmen- 
tation picture. Otherwise, other non-perturbative and 
higher twist effects such as recombination [2(j and the 
influence of radial and longitudinal flow can become im- 
portant and may lead to further modification of the asso- 
ciated yield 2l| . As may be noted from Fig. [2j lowering 
the p?p ssoc ' results in a systematic departure between pre- 
diction and data with increasing centrality. 
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FIG. 2: (Color online) Yields of different flavors of hadrons 
with 6 GeV < Pt soc - < p£ ig ' and 4 GeV < pf soc - < 6GeV, as- 
sociated with a trigger hadron with p^ lg ' > 8 GeV versus the 
centrality of Au+Au collisions at -^s = 200 GeV as compared 
to experimental data [TjJ, see text for details. 

In summary, we have studied the modification of the 
DFFs in both confining and hot deconfined matter due 
to multiple parton scattering and induced gluon radia- 
tion. The modification, follows closely that of SFFs so 
that the associated hadron distributions or the ratios of 
DFFs to SFFs are only slightly suppressed in DIS off nu- 
clei but enhanced in central heavy-ion collisions due to 
trigger bias. With no extra parameters, our calculations 
agree very well with experimental data. Such calcula- 
tions, combined with previous successful comparisons in 
the single inclusive sector [ll| in both cold nuclear and 
hot deconfined matter constitutes a very important test 
of the partonic jet modification formalism and testifies to 



the applicability of perturbativcly calculable jet obscrv- 
ables as probes of the dense matter created in heavy-ion 
collisions. 
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